arXiv: 1111. 1 762vl [astro-ph.SR] 7 Nov 2011 


Mon. Not. R. Astron. Soc. 000 , 1-?? (2002) Printed 9 November 2011 (MN style file v2.2) 


A generalised porosity formalism for isotropic and 
anisotropic effective opacity and its effects on X-ray line 
attenuation in clumped O star winds 

Jon O. Sundqvist 1 *, Stanley P. Owocki 1 , David H. Cohen 2 , 

Maurice A. Leutenegger 3,4 , and Richard H. D. Townsend 5 

1 University of Delaware, Bartol Research Institute, Newark, Delaware 19716, USA 

2 Swarthmore College, Department of Physics and Astronomy, Swarthmore, Pennsylvania 19081, USA 

3 CRESST and X-ray Astrophysics Laboratory NASA/GSFC, Greenbelt, MD 20771, USA 

4 Department of Physics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 
5 University of Wisconsin, Department of Astronomy, Madison, Wisconsin 53706, USA 


Accepted 2011-11-03. Received 2011-11-03; in original form 2011-10-02 


ABSTRACT 

We present a generalised formalism for treating the porosity-associated reduction 
in continuum opacity that occurs when individual clumps in a stochastic medium 
become optically thick. As in previous work, we concentrate on developing bridging 
laws between the limits of optically thin and thick clumps. We consider geometries 
resulting in either isotropic or anisotropic effective opacity, and, in addition to an 
idealised model in which all clumps have the same local overdensity and scale, we 
also treat an ensemble of clumps with optical depths set by Markovian statistics. 
This formalism is then applied to the specific case of bound-free absorption of X- 
rays in hot star winds, a process not directly affected by clumping in the optically 
thin limit. We find that the Markov model gives surprisingly similar results to those 
found previously for the single clump model, suggesting that porous opacity is not very 
sensitive to details of the assumed clump distribution function. Further, an anisotropic 
effective opacity favours escape of X-rays emitted in the tangential direction (the 
Venetian blind’ effect), resulting in a ’bump’ of higher flux close to line centre as 
compared to profiles computed from isotropic porosity models. We demonstrate how 
this characteristic line shape may be used to diagnose the clump geometry, and we 
confirm previous results that for optically thick clumping to significantly influence 
X-ray line profiles, very large porosity lengths, defined as the mean free path between 
clumps, are required. Moreover, we present the first X-ray line profiles computed 
directly from line-driven instability simulations using a 3-D patch method, and find 
that porosity effects from such models also are very small. This further supports the 
view that porosity has, at most, a marginal effect on X-ray line diagnostics in O stars, 
and therefore that these diagnostics do indeed provide a good ‘clumping insensitive’ 
method for deriving O star mass-loss rates. 

Key words: stars: early-type - stars: mass-loss - stars: winds, outflows - radiative 
transfer - line: profiles - X-rays: stars 


1 INTRODUCTION 

Over the past years, it has become clear that in princi- 
ple all standard spectral mass-loss diagnostics of O stars 
are affected by wind clumping, i.e. by the small-scale 
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wind inhomogeneities that should arise naturally from a 
strong, intrinsic instability associated with the radiative 
line-driving of these winds (the line-driven instability, LDI, 
e.g. Owocki et al. 1988). If neglected, such wind clumping 
causes standard diagnostics such as H Q and IR/radio free- 
free emission, which have opacities that depend on the lo- 
cal wind density squared, to overestimate mass-loss rates 
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(for summaries, see Puls et al. 2008; Hamann et al. 2008). 1 
The analysis here examines the degree to which X-ray line 
profiles can provide a mass-loss diagnostic that is relatively 
insensitive to clumping. 

In single O stars without strong magnetic fields, X- 
rays are believed to originate in embedded wind shocks 
associated with the LDI (Feldmeier et al. 1997), and the 
broad emission lines revealed by high-resolution X-ray spec- 
troscopy support this basic scenario (Kahn et al. 2001; 
Cassinelli et al. 2001; Cohen et al. 2006). These X-ray lines 
are often observed as blue-shifted and asymmetric, charac- 
teristics stemming from attenuation by bound-free absorp- 
tion in the bulk wind (MacFarlane et al. 1991). As seen by 
an observer, X-ray photons emitted in the receding part of 
the wind travel farther before escape, and thus undergo more 
absorption, than those emitted in the advancing part. 

For optically thin clumps, the amount of bound-free 
absorption is proportional to the local density, and may 
thereby be used to put additional constraints on mass-loss 
rates. Whilst initial analyses (e.g., Kramer et al. 2003) re- 
quired very low mass-loss rates to reproduce the observed 
X-ray lines, more recent investigations with a better account 
of the wind opacity and line blends (Cohen et al. 2010a, b, 
2011) show that rates inferred from X-ray lines are consis- 
tent with those derived from other diagnostics, if clumping 
is adequately accounted for in the other diagnostics. How- 
ever, a possible shortcoming of these X-ray analyses is the 
assumption that clumps are optically thin, which if not met 
would lead to an overestimate of the wind opacity, due to 
the principal effect of porosity. 

Wind porosity models aiming to calculate X-ray line 
profiles have been developed by, e.g., Feldmeier et al. (2003), 
Oskinova et al. (2004), and Owocki & Cohen (2006) (here- 
after OC06). The first two of these studies assumed the 
clumps to be radially oriented , geometrically thin shell frag- 
ments (‘pancakes’), leading to a distinct anisotropic form of 
the effective opacity. In OC06, on the other hand, the clumps 
were assumed isotropic to impinging radiation. Whereas ge- 
ometrically thin shell structures are indeed seen in one- 
dimensional (1-D) LDI simulations, first attempts to con- 
struct 2-D LDI models suggest that these shells break up via 
Rayleigh- Taylor or thin shell instabilities into small clumps 
of similar angular and radial scales (Dessart & Owocki 2003, 
2005a). But these initial two-dimensional simulations do not 
yet properly treat the lateral radiation transport that might 
couple material, and so the degree of anisotropy of insta- 
bility generated structure in a fully consistent 3-D model is 
still uncertain. 

From the diagnostic side, OC06 (see also Cohen et al. 
2008) argued that for porosity to significantly affect the X- 
ray line formation, unrealistically large porosity lengths, de- 
fined as the mean free path between clumps (see Sect. 3), 
must be invoked. This view is also supported by the above- 
noted recent attempts to derive mass-loss rates directly from 

1 Diagnostic ultra-violet (UV) resonance lines, which have opac- 
ities that depend linearly on density, are directly affected by 
clumping only if individual clumps are optically thick. However, 
recent results indicate that clumps are indeed thick in these 
lines (Prinja &: Massa 2010), which then can lead to reduced line 
strengths and underestimates of mass-loss rates if neglected in 
the analysis (Oskinova et al. 2007; Sundqvist et al. 2011). 
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Figure 1. Schematic of photon of direction n impacting a spher- 
ical clump and a radially oriented (r/|r|) shell fragment. The for- 
mer gives an isotropic effective opacity whereas the latter gives 
an anisotropic effective opacity, since the projected surface area 
depends on direction cosine // = cos 0 = n ■ r/|r|. 

X-ray diagnostics. On the other hand, Oskinova et al. (2006) 
have argued that anisotropic clumps enhance porosity ef- 
fects, and lead to more-symmetric line profiles than if as- 
suming isotropic clumps, in general agreement with X-ray 
observations. 

This paper and its sequel (Leutenegger et al. 2011, 
hereafter Paper II) further examine these issues. Building 
upon previous works by Owocki & Cohen (2001) and OC06, 
Sects. 2 and 3 develop a generalised formalism for synthe- 
sising X-ray lines, including porosity as caused by either 
isotropic (spherical, or randomly oriented) or anisotropic 
(flattened, radially oriented) clumps. In addition, we gener- 
alise our models to treat an ensemble of clumps of some dis- 
tribution in optical depth, rather than retaining the assump- 
tion that all clumps are locally identical. Sect. 4 then system- 
atically examines synthetic X-ray line profiles and analyses 
porosity effects for isotropic and anisotropic clumps, as well 
as for uniform and exponential clump distributions. We dis- 
cuss how the shape of the clumps affects the line profiles 
in cases where porosity is important, and how this may be 
used to put empirical constraints on the wind’s clump geom- 
etry (leaving detailed confrontation with observed spectra 
to Paper II). Sect. 5 presents first X-ray line profiles cal- 
culated directly from LDI simulations, using the 3-D patch 
method first developed by Dessart & Owocki (2002). Sect. 6 
gives a physical interpretation of the analytic porosity mod- 
els presented, showing they can be reconciled with a general 
statistical model derived for stochastic transport in a two- 
component Markovian mixture of immiscible fluids. Finally, 
in Sect. 7 we discuss our results, compare them to other 
studies, and give our conclusions. 


2 OPACITIES IN A CLUMPED HOT STAR 
WIND 

In our phenomenological model, we assume that the opac- 
ities in the bulk wind can be described using a two- 
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Figure 2. Back-lit rendition of randomly generated spherical (upper row) and radially compressed (lower row) clumps in a spherically 
expanding structured wind (columns 2-7), compared to a back-lit smooth wind (column 1) with r* = 1 and the same total mass between 
an onset radius 1.05 R+ and a maximum visualisation radius 5 7?*. The terminal porosity length fooo/f?* increases from left to right, as 
given by the left header, but clump diameters all scale as d c \ = (0.2/y/n)r. The right header gives the total number of clumps. The white 
circle represents the star, which also radiates with the same surface brightness as the background source. Appendix A gives details on 
how these visualisations were generated. 


component medium consisting of overdense ‘clumps’ (de- 
noted with cl) and a rarefied ‘inter-clump medium’. The 
distribution of X-ray emitters in the shock-heated wind is 
described in Sect. 3.4. We neglect the inter-clump medium’s 
contribution to the opacities, an assumption well justified for 
absorption of X-rays in O-star winds, due to the generally 
low X-ray optical depths found for such stars (Cohen et al. 
2010a, 2011). The volume filling fraction of the dense gas is 
fv, thus the local mean density is (p) = fvPch where the 
angle brackets denote spatial averaging. 


2.1 Optically thin clumps 

X-rays emitted in the wind are attenuated by bound-free 
absorption depending linearly on density. The local atomic 
mean volume opacity per unit length is then (\) = k{p), 
with mass absorption coefficient k. By requiring that the 
mean density of the clumped wind be equal to the density 
of a corresponding smooth wind model, i.e. that ( p ) = p sm , 
one immediately recognises the well-known result that for 
atomic processes depending linearly on density, the opaci- 
ties are not directly affected by clumping as long as the in- 
dividual clumps remain optically thin. However, also linear- 
density opacities can via a modified wind ionization balance 
be indirectly affected by optically thin clumping, but mod- 
elling such ionisation equilibria is not a focus of the present 
paper. 

As a comparison, for processes depending on the square 
of the density (e.g. H Q in hot star winds), the opacities are 
always enhanced compared to smooth models, by a factor 
given by the so-called clumping factor f c \ = (p 2 ) / (p) 2 = 
J'y . where the latter equality holds when the inter-clump 
medium is neglected, whereby one obtains 2 (x) oc (p 2 ) = 
fd(p) 2 = /vVL- 


2 In this context, we note that the porosity formalism presented 

in Sect. 3, although developed there for the specific case of X- 
ray line attenuation, is applicable also for continuum processes 
depending on (p 2 ), for example thermal free- free emission, simply 
by exchanging the expression for (y) . 


2.2 Porosity 

It is important to realise that (\) may be accurately used 
in radiative transfer models only in the limit of optically 
thin clumps, r t i < 1. If this condition is not satisfied for the 
investigated process, the radiative transfer becomes more 
complex. For continuum processes such as the attenuation 
of X-rays considered here, optically thick clumps lead to a 
local self-shielding of opacity within the clumps, which in 
turn allows for increased escape of radiation through porous 
channels in between the clumps 3 . The essential effect of 
such porosity is that the ‘effective’ opacity of the medium 
becomes lower than predicted by an optically thin clump 
model. Thus porosity can mimic the symmetrising effects of 
reduced mass loss on the X-ray line profiles. The purpose 
of this paper is to present a formalism for quantifying this 
reduction and to show how clump geometry and distribu- 
tion affect X-ray line profile morphology in the presence of 
porosity. 


3 A POROSITY FORMALISM FOR 
EFFECTIVE OPACITY 

In analogy with the atomic opacity, we may write the effec- 
tive opacity per unit length of a clump ensemble as (e.g., 
Feldmeier et al. 2003) 

Xeff — UclAdP, (1) 

where n c i is the number density of clumps, A c i is the pro- 
jected area (the geometric cross-section) of a clump for di- 
rection n, and P is the probability that a photon impacting 
a clump gets absorbed. This probability obviously depends 

3 For line formation in a rapidly accelerating clumped medium, 
optically thick clumps lead to corresponding velocity gaps, 
through which line photons may escape without ever interacting 
with the material (Owocki 2008; Sundqvist et al. 2010). This is a 
consequence of the Doppler shift, leading to a picture wherein the 
clump length scales can be comparable to (or even larger than) 
the extent of the lines’ resonance zones, which is a limit wherein 
the porosity formalism developed here is not directly applicable. 
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on the clump optical depth, P = P(r c i), which we here 
characterise by the average over all possible chord lengths i 
through the clump for rays of direction n, 

f redAci f f npcidldAci kM c \ , 0<i 

Tcl = A* ~A^i = ^7’ (-) 

where M c \ is the mass of the clump and the last equality 
assumes that k is constant over the clump. 

Another useful quantity is the local mean free path 
between clumps, also known as the porosity length h 
(Owocki et al. 2004), 


h= — 

( 3 ) 

Using this definition, t c i may be written 


Tci = K ^ cl = K(n c \M c \)h = n{p)h = (x)h, 
A-d 

( 4 ) 

whereby 


P(t c i) s Xeff -P(Tcl) 

( 5 ) 

Xeff — , =k , . — 

h (X) Tci 


For P = 1, Eq. 5 returns the atomic opacity independent 
result Xeff = 1/h, demonstrating how the porosity length h 
can also be interpreted as a photon’s mean free path in the 
limit of only optically thick clumps. Further, for a constant 
clump density Eq. 2 yields r c i = Kp c i-£ av , with average chord 
length 4v, which then recovers the commonly used form 
h = £av//v for the porosity length. 


3.1 Isotropy vs. anisotropy 

Note that the porosity length as defined in Eq. 3 is a strictly 
local quantity. And because t c i = (x)h, this means that the 
(an)isotropy of the effective opacity (Eq. 5) depends only on 
A c i, and is independent of the spatial variation of h associ- 
ated with the global wind expansion. Thus, spherical clumps 
(Fig. 1), as well as randomly oriented clumps of arbitrary 
shape, will have an isotropic effective opacity. For these cases 
then, assuming spherical symmetry for the global wind ex- 
pansion, one may set h = h(r) for all directions n impacting 
the clump. 

However, now let us consider a specific wind model 
in which the clumps really are randomly distributed, but 
radially oriented, geometrically thin shell fragments (‘pan- 
cakes’) (Feldmeier et al. 2003; Oskinova et al. 2004). In such 
a model, the projected clump area is A c jn • r/|r| = Ap = i |/u| , 
where we here identify the projected area for a radially di- 
rected ray with the area for an isotropic clump (Fig. 1). 

This implies the effective opacity retains its basic form 
(Eq. 5) also for such anisotropic clumps, but the clump op- 
tical depth becomes larger for oblique rays, 

Td(r,[i) = (x)h{r,p) = (x)h(r)/\p\, (6) 

where the isotropic case is recovered by setting p = 1, 
and the radial dependencies of the mean opacity and di- 
rectional cosine have been suppressed. This paper considers 
only isotropic and ‘pancake’ geometries, but the porosity 
formalism outlined above applies generally to any clump ge- 
ometry described by A c \ = /(n), where / is some function. 

To illustrate such isotropic vs. anisotropic absorbing 
media, Fig. 2 compares a random distribution of spherical 


clumps and radially oriented pancakes, as illuminated by 
a uniform background source (see Appendix A). Note that 
for visual clarity, we extend these visualisations only to an 
outer radius 5 R*. We stress that Fig. 2 is for general il- 
lustration purposes; the uniform background illumination is 
not a distribution of X-ray emitters. Also, we repeat that 
though assumed in the visualisations, spherical clumps are 
actually not a necessity for obtaining isotropic porosity; the 
requirement for this is rather that the clumps be randomly 
oriented, see above. 


3.2 Bridging laws for the effective opacity 


For attenuation with a given local clump optical depth r c \ = 
ro, let us assume the probability of absorption simply takes 
the basic form P = 1 — e~ T ° . While formally exact only in 
cases where all chord lengths across the clump are equal, 
this expression is a suitable approximation that yields with 
Eq. 5 


Xeff _ 1 - e T ° 

(X) ro 


( 7 ) 


This ‘single clump’ bridging law now has the correct values 
in the limiting cases; it returns the atomic mean opacity 
when to <S 1 and is independent of it when to 2S> 1. And as 
discussed in the preceding section, the bridging law equation 
applies for both isotropic and anisotropic effective opacity 
models, however with different expressions for the clump 
optical depth. 

Assuming the effective mean free path scales as Xe® = 
(x) -1 + h, an even simpler ‘inverse’ (or ‘Rosseland’, see 
OC06) bridging law was invoked in OC06, 


Xeff = 1 

(x) 1 + to ’ 


(8) 


which also has the correct optically thin and thick limits. 
By considering the direction dependent mean free path, this 
inverse bridging law can be realised also for anisotropic mod- 
els. 

In OC06, the practical motivation for invoking Eq. 8 
was because the optical depth integral for X-ray attenuation 
could then be solved analytically. However, Sect. 6 shows 
it also happens to represent the bridging law that follows 
from assuming the local clump optical depth distribution 
function obeys Markovian statistics. Thus, the single-clump 
and inverse bridging laws (Eq. 7 and Eq. 8) differ in that the 
former assumes all clumps have the same local optical depth 
(for a given direction), whereas the latter averages over an 
exponential distribution in t c i. 


3.3 Velocity stretch porosity 

OC06 assumed that the porosity length scales with the lo- 
cal radius, but for mass-conserving clumps, such a scaling 
is only appropriate for isotropic expansion. For clumps re- 
leased into a radially expanding stellar wind, the wind ac- 
celeration will ‘stretch’ the clump spacing in proportion to 
the wind velocity (Feldmeier et al. 2003). The analysis here 
assumes this velocity stretch form for both isotropic and 
anisotropic porosity. 

This distinction is most easily seen for the radially frag- 
mented shell model, in which the average radial separation 
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between two shells is A r = h(r). For shells moving radi- 
ally according to a ‘/3-velocity law’, w(r ) = v(r)/voo = 
(1 — R+/r)P, where 7?* is the stellar radius, the separa- 
tion is h(r) = /ioow(r), with the parameter hoo representing 
the asymptotic radial separation as w — > 1. For simplic- 
ity, this paper assumes the prototypical value /3 = 1. The 
quantity hoo/voo represents the average time between two 
consecutive shell passings at a fixed radial point in the wind 
(Sundqvist et al. 2010), which may also be interpreted as the 
inverse of a ‘fragmentation frequency’ no (Oskinova et al. 
2004, 2006). For example, hoo = R * gives a fragmentation 
frequency no = Voo/R* that is equal to the inverse of the 
wind flow time. 


3.4 X-ray line transfer in porosity models 

To compute X-ray emission line profiles, we solve the stan- 
dard formal integral of radiative transfer, using a customary 
( p , z) coordinate system and following the basic procedure 
described in Owocki & Cohen (2001) for the distribution of 
X-ray emitters. Since our primary interest here is the shapes 
of the lines, all resulting flux profiles have been normalised 
to a unit maximum, jrnorm _ jr^yjyjax^x). For simplicity 
we assume that the X-ray emission begins at a certain onset 
radius Ro = 1.5 7?*, and is constant beyond it. This onset 
radius is consistent with that typically predicted by con- 
servative, self-excited LDI simulations (Runacres & Owocki 
2002), and both Ro « 1.57?.* and a constant X-ray Hlling fac- 
tor are supported by observations (Leutenegger et al. 2006; 
Cohen et al. 2006, 2011). Moreover, both LDI simulations 
and the observed X-ray luminosities indicate that only a 
very small mass fraction, less than 1 %, of the stellar wind is 
shock-heated to X-ray emitting temperatures at any given 
time. Simulations and observed lack of X-ray variability fur- 
ther indicate that there are numerous sites of X-ray emis- 
sion distributed throughout the wind, justifying our assump- 
tion of a smoothly distributed X-ray emitting plasma above 
the onset radius. Note though, that Ro does not necessar- 
ily equal the clump onset radius 7? c i, which observations 
typically indicate is located much closer to the photosphere 
(Puls et al. 2006; Cohen et al. 2011, see also Fig. 2 for a vi- 
sualisation). But in the velocity stretch porosity models, we 
have verified that the exact value of this R c \ is not important 
for the resulting X-ray line profiles. 

The absorption of X-rays emitted at position Z e along 
a ray with impact parameter p is given by the optical depth 
integral 


r(p,2e)= / Xes{z,p)dz, (9) 

J Z e 

where the effective opacity accounts for any porosity. The 
opacity in a smooth or optically thin clump model, due 
purely to the atomic mean opacity (x), is proportional to 
the mass-loss rate M of the star and here characterised by 
a fiducial optical depth r* = Mk/IAttR+Voo), with wind ter- 
minal speed Voo. To evaluate Eq. 9 for our porosity models, 
the only additional input parameter required is the poros- 
ity length h (see Eqs. 4-8). This holds for isotropic as well 


as anisotropic effective opacity models, and for the single- 
clump as well as the inverse bridging laws. 4 


4 X-RAY LINE PROFILES FROM ANALYTIC 
POROSITY MODELS 


Fig. 3 displays synthetic X-ray line profiles calculated using 
the four possible combinations of isotropic vs. anisotropic 
effective opacity and single-clump vs. inverse bridging laws. 
The figure clearly shows that profiles calculated using the 
two different bridging laws are very similar (see also Fig. 5, 
as well as Fig. 1 in OC06), despite representing two very 
different clump optical depth distributions (Sect. 3.2). This 
indicates that the effects of porosity on X-ray line profiles 
are not very sensitive to the specific local distribution in r c i. 
We discuss this important result further in Sect. 6. 

The second key feature of Figs. 3 and 5 is the promi- 
nent ‘bump’ visible close to line centre in profiles calculated 
with anisotropic effective opacity. Conceptually, we may un- 
derstand this as a Venetian blind’ effect (Fig. 4, see also 
Feldmeier et al. 2003); since the fragmented shells are radi- 
ally oriented, the blinds are closed for radial photons, but 
open up for more tangential ones. This leads to increased 
escape for photons emitted close to line centre, since the 
line emission wavelength scales with direction cosine p as 
x = — pw (e.g., Owocki & Cohen 2001). 

Another way to look at this effect is to consider the 
optical depth integral for anisotropic effective opacity in the 
Td » 1 limit, 


r(p,z) 



Xeff(2, 




(10) 


which shows that, since dr = pdz, the optical depth in this 
limit is set simply by counting up the number of porosity 
lengths. In the plane-parallel limit of radially oriented, geo- 
metrically thin but optically thick fragments, all tangential 
(p = 0) photons would escape. However, due to sphericity ef- 
fects (i.e. , that p increases as the photon propagates through 
the wind) also photons emitted initially in the tangential 
direction will suffer some absorption (see Fig. 4). Thus the 
end result is not complete transmission, but a characteristic 
bump stemming from the reduced integrated optical depth 
for photons emitted around x ~ 0. For isotropic porosity, on 
the other hand, no p factor enters in Eq. 10, and therefore 
no bump appears in these profiles. 

This quite distinct and systematic difference in the 
shape around line centre between models with isotropic and 
anisotropic effective opacity is a key result of the present 
analysis. Indeed, one can us this difference to set empirical 
constraints on the clump geometry by confronting synthetic 
X-ray spectra with observed ones, as will be done in Paper II. 

Generally, Fig. 3 confirms earlier results by OC06 that 
in order to achieve a significant effect on the profiles, rather 
large porosity lengths, hoo > 7?*, are required. However, 


4 Source codes to all X-ray line porosity models presented 
in this section are publicly available at (package windprof) 
http: / /heasarc. gsfc.nasa.gov/docs/xanadu/xspec/models/wind 
prof.html, where models for the broadband absorption of X-rays 
(package WINDTABS, Leutenegger et al. 2010) also can be found. 
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Isotropic opacity Anisotropic opacity 


hJR.=0.0 


hJR.=0.5 


hJR.=1 .0 


hJR.=2.0 


hJR.=5.0 


Xeff/ < X >= (1-e" T )/x 




Xef/ < 5C >= 1/(1+x) 


Xeff/ < X > =( / >-e' t )/x 



Figure 3. Synthetic X-ray line profiles for the ‘single clump’, Xeff / (x) = (1 — e — r )/r, and the ‘inverse’, Xeff/(x) = 1/(1 + T ); effective 
opacity bridging laws and for different porosity length parameters hoo/R *, as labelled. All models assume an onset radius for the X-ray 
emission Rq = 1.5/J*. The abscissae display the dimensionless wavelength x = (A/Ao — l)c/uoo, with Ao the line-centre wavelength, and 
the ordinates display the normalised flux. Black, red, blue, and green dashed lines have an optical depth parameter r* = 0.01, 1,5, 10, 
respectively; for a non-colour separation, an increased T r means a more blue-shifted peak flux. Note that we have set hoo = 0 (the 
uppermost panel) to be equivalent to assuming only optically thin clumps. 


Fig. 5 reveals that for a value of t* = 2.5, representa- 
tive of the prototypical O supergiant £ Pup (Cohen et al. 
2010a), the anisotropic porosity model displays significantly 
higher flux around line centre than the other profiles, also 
for /loo = f?*. Further comparisons show that, for the pa- 
rameters used in Fig. 5, the isotropic porosity model is well 
matched by a corresponding optically thin clump model with 
r* reduced by ~ 20 percent (a detailed quantification of this 
trade-off will be provided in Paper II). In contrast, such a 
simple optical depth reduction does not reproduce the line- 
centre region in anisotropic porosity models, i.e. there is 
no simple trade-off, or degeneracy, between mass-loss rate 
and anisotropic porosity. Thus, whereas it will be difficult 
to distinguish between optically thin clumps and moderate 
isotropic porosity with somewhat higher r*’s, careful line- 
fitting to observations should be able to identify, or refute, 
an anisotropic porosity. 


in a large portion of the wind, the profile shapes do indeed 
become quite independent of atomic opacity. In the bottom 
panel of Fig. 3, profiles computed using an anisotropic ef- 
fective opacity are actually very similar to the flat-topped 
profiles stemming from computations without any X-ray ab- 
sorption, whereas profiles computed using an isotropic ef- 
fective opacity retain a certain degree of asymmetry. This 
result differs from that found in OC06, wherein very near 
symmetry was achieved for isotropic profiles with very large 
porosity lengths, and comes about due to the h oc v scaling 
adopted here, which implies significantly shorter porosity 
lengths in the lower wind regions than the h oc r used in 
OC06. However, for more moderate values of /loo (top three 
panels in Fig. 3), all profiles are opacity dependent, as ex- 
pected since in the formation of these lines, porosity is a 
secondary effect. 


Finally, for very large porosity lengths, giving r c i 1 
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Figure 4. Illustration of the ‘Venetian blind’ effect seen in poros- 
ity models using an anisotropic effective opacity. The dashed ar- 
rowed lines represent two different p-rays and the observer is as- 
sumed to be located at Zoo- 



Figure 5. Line profiles for hoo = R* and r* = 2.5, using different 
effective opacity laws, as labelled. 



Figure 6. Line profiles for r* = 2.5, calculated from a smooth 
CAK model and structured LDI models with patch sizes 1 and 3 
degrees (see text), as labelled. 


5 X-RAY LINE PROFILES FROM LDI 

SIMULATIONS 

The above simplified analytic models demonstrate that the 
effect of optically thick clumps on wind continuum absorp- 
tion can be well characterised in terms of the wind porosity 
length, parametrised here by the asymptotic value /too and 
assuming a radial variation set by the velocity-stretch form, 
h(r) = hoo(l — Rt/r ) (Sect. 3.3). 

But inspection of LDI simulations suggests a substan- 
tially steeper radial variation for the separation between 
instability-generated wind clumps. Even in 1-D models in 
which the separation can become of the order of a stellar 
radius in the outer wind (implying hoo « R„), the initial 
clump structure formed near the onset radius r « 1.577, 
tends to have a much smaller separation scale, of the order 
of the Sobolev length L so b « (vth/v)R* « 0.0177*, where 
Vth is the ion thermal velocity. The sharp increase in separa- 
tion from this onset comes not just from velocity stretching 
from the overall wind acceleration, but also from collisional 
merging of clumps with substantial radial velocity disper- 
sion. For a given asymptotic porosity length hoo, LDI mod- 
els thus tend to have smaller inner-wind porosity lengths 
than assumed in the simple velocity stretch scaling. Since 
it is in this inner wind region that clumps can become op- 
tically thick for X-rays, this suggests that LDI models will 
show even weaker porosity effects than implied by the an- 
alytic profiles shown in Figs. 3 and 5. To demonstrate this 
explicitly, we now present some first sample calculations of 
X-ray profiles computed from 1-D LDI simulations that are 
phased randomly among 3-D patches of a parametrised an- 
gular size (as first developed in Dessart & Owocki 2002). 
The details of this patch geometry as implemented in our 
radiative transfer code are given in Sundqvist et al. (2011). 
Here we adapt this code to synthesise X-ray line profiles by 
making the following assumptions: The X-ray emission is 
assumed to have a fixed spatial form independent of wind 
structure, scaled in proportion to the density squared of a 
smooth CAK wind and with an onset radius Ro = 1.5 77*. 
The bound-free absorption is then calculated directly from 
the structured LDI simulation presented in Sundqvist et al. 
(2011) (computed following Feldmeier et al. 1997). 

Fig. 6 plots X-ray line profiles for the same wind optical 
depth used for analytic models in Fig. 5, namely r* = 2.5. 
The curves compare a smooth CAK model to structured 
LDI models with patch sizes of 1 and 3 degrees. The overall 
shapes agree well with corresponding non-porous analytic 
models, except for small differences due to the fact that the 
CAK velocity law does not exactly follow the phenomeno- 
logical /3 = 1 law. But the key point is that the LDI profiles 
are very similar to the CAK profile, implying little porosity 
effect; however, we note that for a 3-degree patch size, there 
is a small, but noticeable, ‘bump’ around line centre, pre- 
sumably associated with the anisotropic (pancake) nature of 
the clumped structure in such models. 

These results confirm that porosity is as marginal 
in LDI simulations as in corresponding analytic stretch- 
porosity models with comparable asymptotic clump sep- 
arations (i.e., hoo ~ R *)■ Moreover, in 2-D LDI models 
(Dessart & Owocki 2003, 2005a), clumps can also be bro- 
ken up by shearing and associated effects, leading generally 
to more, smaller, and less optically thick clumps, charac- 
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terised by even smaller porosity lengths. Future work will 
implement our radiative transfer tools also in such genuinely 
multi-D instability models. 


6 A PHYSICAL INTERPRETATION OF THE 
EFFECTIVE OPACITY BRIDGING LAWS 

As already noted, a basic difference between the two bridg- 
ing laws adopted in this paper is that one (Eq. 7, the ‘single 
clump’ law) assumes a locally constant clump optical depth 
for all clumps, whereas the other (Eq. 8, the ‘inverse’ law) 
represents a certain distribution in r c i. Before discussing 
porosity effects on the synthetic X-ray line profiles in Sect. 7, 
this section examines the nature of this distribution. 


6.1 Distribution laws for r c i 

To investigate how a distribution of clump optical depths 
affects effective opacity scalings, let us assume that the ratio 
of the effective to the mean opacity scales as the ratio of an 
effective clump optical thickness to its mean, 

Xeff Teff / 0 °° T f( T ) hZ V~ dT nn 

<X> <r> / 0 °°r/(r)dr ’ 1 j 

where we for clarity have dropped the indices on the clump 
optical depths. Eq. 11 introduces /(r), the normalised 
distribution function of clumps, and r e fr, the distribution 
weighted mean of the clump optical depth. Selecting a 
weighing function (1 — e _T )/r ensures that the single clump 
bridging law (Eq. 7) is recovered from Eq. 11 when the dis- 
tribution function is a Dirac delta function, S(t — to). 

Let us now choose a specific distribution function of the 
exponential form 
p -t/t 0 

fir) = , (12) 

TO 

which has a mean value to. Using Eq. 12 in Eq. 11 yields 
directly the inverse bridging law Eq. 8, but with to now only 
representing the mean clump optical depth, rather than a 
unique one as in the exponential bridging law Eq. 7. Again, 
despite the large difference between this clump distribution 
and the one assuming a constant t, the two bridging laws 
give similar results, as demonstrated in Sect. 4. 


6.2 Connecting Markovian statistics to 
exponentially distributed clumps 

We now show that the bridging law resulting from this ex- 
ponential distribution turns out to be a special case of a 
general scaling derived for a stochastic mixture of two fluids 
that follow Markovian statistics (e.g., Levermore et al. 1988; 
Pomraning 1991). The Markov assumption is that the future 
state of the system only depends on its present state, and 
not on its history. Along any given ray through the medium, 
if the fluid is of component 0 at location s, the probability 
of it being of component 1 at s + ds is Po.ids, where Po,i 
is independent of how far back along the ray the last tran- 
sition (from fluid 1 to 0) occurred. Under this assumption, 
the length scales traveled within the fluid components are 
random variables described by Poisson distributions, with 


Pop identified as the inverse of £o, the mean distance a pho- 
ton travels along the ray in fluid 0 before finding itself in 
component 1. A similar definition applies for l\. 

For such a two-component Markov model with opacities 
Xo and xii Levermore et al. (1988) derived for the effective 
opacity in the pure absorption case 


Xeff = 


(X) + XoXilc 
1 + (poXi +PiXo)£c’ 


(13) 


where pi = ti/{l o + £i), (x) = PoXo + XiPi is the mean 
opacity, and i c = £o£i/(lo + £i) is the correlation length 
(Pomraning 1991). Identifying the clumps in our model with 
component 1, and assuming the inter-clump medium to be 
void (xo = 0), we find for this ‘clump-|- void Markov model’ 

1 1 1 (») 


Xeff _ 

(x> l -poxi4 


1 + (x)Pof'o l + (x)h' 


Here po represents the probability that a photon is in the 
void medium, while la is the distance the photon travels 
in the void before encountering a clump; the product polo 
thus represents the photon mean free path in the case of 
optically thick clumps, which is also the porosity length h, as 
given by the final equality in Eq. 14. Comparison with Eq. 8 
then shows that the effective opacity bridging law for an 
exponential clump optical depth distribution is equivalent 
to that for this statistical clump+void Markov model. 

Indeed, recalling that the Markov transport is defined 
along a given ray, we may make the same identification for 
the anisotropic porosity model, with the photon mean free 
path along the ray then being scaled by l/|/r|. 


6.3 Exponentially truncated power-law 
distributions 

While there are not many observational constraints on 
the distribution of clumps in a hot star wind (see how- 
ever Lepine & Moffat 1999; Dessart & Owocki 2005b), the 
above identification with the Markov model at least places 
our porosity models on a robust and well-known statisti- 
cal ground. The inverse bridging law should therefore be an 
appropriate standard choice for porosity applications such 
as the X-ray line formation considered here, but perhaps 
also for, e.g., porosity moderated continuum driven wind 
models of stars formally exceeding the Eddington lumi- 
nosity, as investigated by Owocki et al. (2004). Indeed, al- 
though not explicitly studied in that paper, we note that 
the Markov model represents a special case of the exponen- 
tially truncated power-law distribution of clumps considered 
in Owocki et al. (2004), namely the one with power index 
a p = 2. Thus, reasonable extensions of the two canonical dis- 
tributions studied here could readily be done by using some 
other power index variant given in Owocki et al. (2004). 


7 DISCUSSION AND CONCLUSIONS 
7.1 Isotropy or anisotropy 

Let us next compare our analysis to that by Oskinova et al. 
(2006). These authors also pointed out the differences be- 
tween isotropic and anisotropic effective opacity and car- 
ried out a comparison, however only for the specific case 
of t* = 10 and a fragmentation frequency (see Sect. 3.3) 
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no = 1.4 x 10~ 4 s^ 1 . Taking the parameters for the 0 su- 
pergiant C Pup adopted in Oskinova et ah, this corresponds 
to /loo = Vao/no = 4.247?*, which in turn roughly corre- 
sponds to the bottom panel in our Fig. 3. Indeed, the pro- 
files displayed in that panel agree well with those in Fig. 16 
of Oskinova et al.; both figures illustrate that for such very 
large porosity lengths, profiles computed using anisotropic 
effective opacity are nearly symmetric (Sect. 4). 

This comparison suggests an overall good agreement 
among the results found by the different groups. But as 
shown in Sect. 4, anisotropic porosity line profiles, with 
their characteristic ‘bump’ at line centre, are qualitatively 
different than isotropic porosity or optically thin clump- 
ing profiles. Thus the good statistical fits presented for £ 
Pup by Cohen et al. (2010a), without invoking porosity, 
seem somewhat contradictory to the good visual fits pre- 
sented by Oskinova et al. (2006), using models with mod- 
erate anisotropic porosity hoo ~ R+- Paper II will further 
examine and quantify these differences between anisotropic 
porosity on the one hand, and isotropic porosity or optically 
thin clumping on the other. 


7.2 Is porosity important for X-ray line mass-loss 
diagnostics? 

The LDI simulations presented in Sect. 5 indicate small 
porosity lengths and negligible porosity effects on X-ray line 
profiles. Such small porosity lengths also have some indirect 
empirical support. Namely, the mass-loss rate derived for £ 
Pup by Cohen et al. (2010a), without invoking porosity, is 
only marginally lower than the upper limit mass-loss rate 
derived by Puls et al. (2006), by assuming an unclumped 
outermost radio emitting wind, while allowing for clumping 
in the intermediate and lower wind. Because of the trade-off 
between porosity and mass-loss rate then, if porosity lengths 
large enough to significantly affect the X-ray line profiles 
were to be adopted, /loo > 7?*, the inferred X-ray mass-loss 
rate would be higher than this upper limit. That is, such 
multiwavelength considerations indicate that a significant 
porosity effect on X-ray based mass-loss rates is incompati- 
ble with diagnostic results from other wavebands. 

Overall, we thus conclude that porosity effects on X- 
ray line profiles are likely to be, at most, a marginal effect 
in typical O stars. This is supported also by the low opti- 
cal depths found for £ Pup as well as for the even denser 
wind of HD93129A (Cohen et al. 2010a, 2011). Since most 
O stars will have characteristic r*’s significantly lower than 
these, porosity effects should be negligible. The upshot is 
that X-ray line analysis may indeed provide the best avail- 
able ‘clumping insensitive’ diagnostic of O star mass-loss 
rates. 


ACKNOWLEDGMENTS 

This work was supported in part by NASA ATP grant 
NNX11AC40G. D.H.C. acknowledges support from NASA 
ADAP grant NNX11AD26G to Swarthmore College. We 
thank A. Feldmeier for providing the instability simulations 
discussed in Sect. 5, and for suggesting the exploration of 
Markov models discussed in Sect. 6. 


APPENDIX A: GENERATION METHOD FOR 
CLUMPED MEDIUM ILLUSTRATIONS 


The clumped medium illustrations in Fig. 2 were generated 
by following the radial expansion of mass-conserving clumps. 
As with the analogous illustration in Fig. 3 of OC06, we 
assume the clump scale l increases in proportion to the local 
radius, 

l ( r ) = i*jir ■ ( A1 ) 


However, instead of the OC06 assumption of a purely 
isotropic (‘Hubble-law’) velocity expansion v ~ r, we now 
use a standard (3 = 1 wind velocity law. For clumps of 
projected area A c \ = l 2 oc r 2 and local volume density 
n c i oc 1 /(vr 2 ), this gives the associated radial variation of 
the porosity the desired ‘velocity-stretch’ form, 


h(r) = 


ri ci A c i 


= h c 


v(r) 


= hoo( 1 — R*/r) . 


For specified clump parameters Z* and hoo, 
ber density is thus given by 


n c i(r ) = 


VaoRZ 


hA c \ hooll r 2 v(r) ' 


(A2) 

the clump num- 
(A3) 


Note that, unlike the OC06 isotropic expansion model, the 
clump volume filling factor in this velocity-stretch scaling 
is not constant, but varies spatially as /v oc n c i/ 3 oc r/v(r), 
which actually is quite consistent with derived observational 
constraints (e.g., Puls et al. 2006). 

The cumulative number of clumps up to a radius r 
above the clump onset radius 7? c i is 


N(r) 


nr 

4-7T / n c ir' 2 dr' 

Jr ci 

47r7?* r r — 7?d / r — R* V 

h x ll V Rt + \R c i- R„) . ' 


(A4) 


For a specified outer radius 7? max , the total number of 
clumps is TVtot = TV(7? max ). Since Fig. 2 uses fixed param- 
eters U/R* — 0.1, R c i/Rt = 1.05, and 7? max /7?» = 5, this 
number scales with l//ioo, and ranges from TVtot = 1308 
for the largest porosity length hoo/R * = 8 in the right- 
most column, to Ntot = 41 881 for the least porous case 
hoo/R * = 0.25 in column 2. 

A random set of TVtot clumps with the required statisti- 
cal distribution in radius can now be drawn simply by gen- 
erating a set of pseudo-random numbers TZi over the range 
[0, 1], and inverting the normalised cumulative distribution 
function F(r) = N(r)/Ntot = Ri to find the radius r, 


r(TZi) = R, (1 + ProductLog[exp(C c i + C max 7?.i)]) , (A5) 

with the constants given by 

G c i — R c i/R, — 1 T ln(7? c i/7?* — 1) (A6) 

C max = (7? max - R c \)/R* + In ( ^ ax ) • (A7) 

V ltd — lt» / 

Likewise, we use the assumed statistical spherical symmetry 
and generate the clump angle coordinates in azimuth cj> and 
colatitudinal cosine p through additional pseudo-random 
numbers, 


(pi = 2nTZ i+Ntot ; pi = 21Zi + 2N tot — 1 • (A8) 
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In this system it is most convenient to assume the clumps 
are viewed from above the ^ — 1 pole, with thus fi > 0 
(/r < 0) representing clumps in the foreground (background) 
hemisphere. 

With the random set of clump positions in hand, the 
clump sizes are set by Eq. (Al), with all cases in Fig. 2 
using U/Rt = 0.1. To give the associated projected area 
l 2 = A c i = 7rdj?i/4, the clump diameters are set to d c i = 

2 Z/v^F. 

For the spherical clumps in the upper row of Fig. 2, the 
transparency of each individual projected clump disk area 
is set by exp(— t c i), where the surface-averaged clump opti- 
cal depth is T c i = TthaoRt/r 2 . For the radially compressed 
clumps in the lower row, the associated clump optical depths 
are increased by l/|/z| ; their projected areas are reduced 
through foreshortening their radial extent by a factor 
while keeping their radially perpendicular extent equal to 
the local clump diameter d c \. 

Finally, clumps in the back hemisphere (fi < 0) that 
are behind the star (with T\J 1 — fj? < R,) are simply not 
drawn. This effectively means the clumps directly in front of 
the star appear as if illuminated by a stellar surface bright- 
ness equal to the back illumination of the clumps outside 
the stellar limb. 
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